Sensitivity of measured fission yields on prompt-neutron corrections 
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The amount of emitted prompt neutrons from the fission fragments increases as a function of 
excitation energy. Yet it is not fuUy understood whether the increase in v {A) as a function of E-n 
is mass dependent. The share of excitation energies among the fragments is stiU under debate, 
but there are reasons to beheve that the excess in neutron emission originates only from the heavy 
fragments, leaving vught [A) almost unchanged. In this work we investigated the consequences of a 
mass-dependent increase in h' [A) on the final mass and energy distributions. The assumptions on 
V {A) are essential when analysing measurements based on the 2E-technique. This choice showed 
to be significant on the measured observables. For example, the post-neutron emission mass yield 
distribution revealed changes up to 10-30 %. The outcome of this work pinpoint the urgent need 
to determine v {A) experimentally, and in particular, how v [A) changes as a function of incident- 
neutron energy. Until then, many fission yields in the data libraries could be largely affected, since 
they were analysed based on another assumption on the neutron emission. 



I. INTRODUCTION 

Accurate fission-fragment (FF) yield measurements are 
essential both for reactor- and model calculations. They 
are needed in various stages of the nuclear fuel cycle, e.g. 
reactor criticality- and safety calculations, burnup deter- 
mination, decay-heat calculations and processing of spent 
nuclear fuel [T]. Fission yields are also important for the 
development of fission models. Countless yield measure- 
ments have been performed utilizing different techniques. 
This study focuses on the widely used double energy tech- 
nique (2E) in which the mass yields are determined by 
measuring the FF kinetic energies. This technique suf- 
fers from the vital need to assume the neutron emission 
on an event-by-event basis. In this work we show that 
this assumption has a rather strong implication on the 
measured yields. 

The accelerated FF remove the major part of their 
excitation energies by emitting prompt neutrons. The 
prompt-neutron multipHcity v depends on the fragment 
mass, kinetic energy and excitation energy in the fission- 
ing system. The mass dependence is governed by the 
deformation of the fragments and by shell effects. Due 
to the doubly magic nucleus at ^ = 132, the neutron 
emission is highly suppressed. The resulting v {A) shape 
has a "saw-tooth" like structure and is shown for in- 
stance in Figs. l[a,b) [HE]. Measurements of the to- 
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tal average neutron multiplicity (ftot)) e.g. Ref. [1] for 
^■^■*U (n, /), showed an increasing trend as a function of 
£;„. At ^n ~ MeV, v{A) is roughly 2.4 and it be- 
comes roughly 3.0 at i?n = 5 MeV. In various models 
and experiments, these extra neutrons are assumed to 
be independent on the fragment mass [5Hll). Therefore, 
V {A) was increased almost equally for the light and heavy 
fragments. However some measurements found that the 
extra neutrons are predominantly emitted from the heavy 
fragments, e.g. in ^ar^p („ j) g] and ^ssu (n, /) 
(see Figs. l^a,b)). These experiments used the (2E-2v) 
technique which is not confined by the neutron emis- 
sion assumption since the FF velocities are also mea- 
sured. Both pre- and post neutron-emission distribu- 
tions were determined, which allowed for an estimation 
of v{A). Recently, theoretical models explained this in- 
crease in i^hoavy {A) based on the energy-sorting mecha- 
nism [T^. In the pre-scission state, the fragments are still 
connected through the neck and can exchange excitation 
energy. The temperature is defined as the mean excita- 
tion energy per degree-of-freedom, and depends on the 
fragment mass as T oc A~'^/^ [T^. The light fragments 
have higher temperatures than the heavy ones. Because 
of the thermal contact, the energy will flow to the heavy 
fragments. Due to pairing correlations the nucleons rear- 
range so that the temperature stays constant during the 
entire process, despite the added excitation energy to the 
heavy fragment. The energy sorting is believed to behave 
similarly for neighbouring nuclei. The GEF code uses the 
energy sorting and can be used to estimate the expected 
behaviour I4j. A GEF calculation was performed for 
^■^^U (n, /) and can be seen in Fig. Ijc). 
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II. THE 2E-TECHNIQUE 

The present work aims to study the sensitivity of the 
measured yield distributions on the assumption of v {A) . 
Experimental data on ^'^'^U (n, /) were used at i?n = 4 and 
5 MeV [7j . It is still argued whether the change of v (A) 
at different E^ is mass-dependent. However if the neu- 
tron emission proves to increase only from the heavy frag- 
ments, the measured yields based on the 2E-technique, 
might be affected if they have been calculated using the 
average increase model. In the 2E-technique both FF ki- 
netic energies are measured in combination usually with 
their emission angles relative to the incoming beam. Con- 
servation of momentum and mass are used to iteratively 
calculate the pre-neutron emission masses (for details on 
the kinematics see Ref. [7]). The detected post- neutron 
energies [E\,o^^ are converted to pre-neutron energies 
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In this approximation the recoil from the neutron emis- 
sion to the accelerated fragments is neglected. In most 
cases, the neutron emission is not measured so an approx- 
imation is needed. On average, the neutrons are isotrop- 
ically emitted in the CM system. Moreover the added 
recoil is very small due to the mass difference of the neu- 
tron and the fragment. Unless measuring the neutron 
emission on an event-by-event basis, these assumptions 
are unavoidable even when using the 2E-2v method. The 
iJp^^ are converted into the centre-of-mass system and 
used to calculate the pre-neutron emission masses: 



Al,2=AcNi?2Tx(i?r+i?2™)" 



(2) 



were (1,2) denote the two FF. The neutron multiplicity is 
updated in each iteration based on the newly determined 
masses and energies. The neutron emission has to be 
parametrized as a function of mass and TKE. Moreover, 
the total neutron emission has to correspond to i^tot from 
Ref. [1] . Two approaches (see Fig. 2 ) of the correction 
were tested individually and compared side-by-side: 

Average method (AV) : v {A) was assumed to in- 
crease for all masses. The original v (A) was multi- 
plied by a factor (a\) and adjusted (weighted with 
the mass yields) to correspond to the expected j>tot • 

Heavy method (HE): v {A) increases for A > 120. 
The original sawtooth curve for the heavy masses 
was multiplied with a factor (02) to give the to- 
tal expected neutron emission. For A < 120, the 
distribution was kept unchanged at i^th (A) . 

When applying the AV method e.g. on the 5 MeV case, 
both the light and heavy fragments get an increase of 
roughly 0.3 neutrons. Based on the HE method instead, 
the heavy fragments emit all extra 0.6 neutrons. In both 
methods, a is fine-tuned to give Ptot ~ 3 neutrons. 
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FIG. 1. (a) iy{A) for '^'^Np (nj) at different En [1]. (b) 
u{A) for ^^^U(n, /) at different En [3]. (c) GEF calculation 
performed for ^^''U (n, /) at En = and 5 MeV. 
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FIG. 2. !?234 (A) for En = and 5 MeV deduced from the 
average of 233,235-[j jjg rpj^g original 1/234 (A) shape was mul- 
tiplied by a factor (a) and weighted with the mass distribution 
to give i/tot at _En = 5 MeV. Two different assumptions were 
compared for the increase in i? as a function of En- The first 
using Qfi acting over all masses, the second using 02 acting 
on the heavy fragments only (A > 120). In order to yield the 
same i/tot at En = 5 MeV, 02 > ai. 
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FIG. 3. (a) The post-neutron mass distribution for ■^^*U (n, /) at En = 5 MeV, for both applied corrections (AV and HE), (b) 
The ratio between the distributions, (c) The difference in (TKE) as a function of pre-neutron emission mass, (d) The difference 
in (-Epost) as a function of post- neutron emission mass. Only statistical errors are shown. 



III. IMPACT OF THE DIFFERENT 
ASSUMPTIONS 

The difference in neutron emission per fragment had 
a clearly visible impact on the distributions. The two 
mass peaks moved closer to each other by using the HE- 
method instead of the AV-method. This was observed 
both in the pre-neutron emission and the post-neutron 
emission distributions. However the effects were at least 
twice as pronounced in the latter case. The impact is in- 
directly affecting the pre-neutron masses during the itera- 
tion calculations. But it directly affects the post-neutron 
masses since those are determined by a direct subtrac- 
tion Apost = ^prc — i^- Fig. 3 shows some examples of 
the effects observed. In (a) the post-neutron mass distri- 
butions for En = 5 MeV are plotted. The distributions 
show an overall shift and a small change in the heavy 
peak width. Specific masses experience large differences, 
especially around A w 132. The ratio between the two 
distributions is plotted in Fig. 3jb). As seen, the differ- 
ences are about 10-15 % for masses with 0.04 absolute 



yield. For masses with 0.03 absolute yield, the differ- 
ences reach between 20-30 %. These large discrepancies 
will have a considerable influence on reactor- and model 
calculations. The accuracy often requested is far below 
these levels. Furthermore, it was found that all effects in- 
crease with the incident-neutron energy therefore larger 
corrections are needed when using yields at higher excita- 
tion energies. The energy distribution is also affected as 
seen in Fig. 3jc). The change in TKE distribution reach 
up to ATKE == 0.5 MeV for masses around A ^ 132. 
For higher asymmetries, the differences are compensated 
between the AV and HE distributions. The single frag- 
ment kinetic energies are plotted in Fig. 3(d) as a func- 
tion of Apost- The heavy fragments overall have larger 
kinetic energy when using the AV-method. The light 
fragments change however differently: The HE-method 
gives higher kinetic energies but most notably at higher 
asymmetry. The observed changes are obviously mass- 
dependent which increases the complexity in the needed 
corrections. For the full details on the observed differ- 
ences, see Ref. [TBI. 
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The average changes in TKE and mass distribution 
were studied. The TKE obtained from the HE-method is 
smaller than the AV-method and can be estimated by: 



TKEhe ~ TKEav - 0.038 x £:„ (MeV) (3) 

The pre-neutron mass distribution changes roughly by: 

«°>HE « «°>Av - 0.065 X S„ (u) (4) 

The post-neutron mass distributions showed double the 
effect compared to the pre-neutron case. At E^ — 5 MeV, 
the change in average mass is about 0.7 amu: 



(Ar*)HE « (^r'>AV - 0.135 X E^ (u) 



(5) 



Eqs. (3, 4, 5) may be used as a first correction to the 
expected shift. However, average shifts are probably not 
sufficient to account for the total changes, since the dif- 
ferences seen in Fig. 3(a-d) are mass-dependent. 



IV. CONCLUSIONS 

In this work, we investigated the impact of differ- 
ent D (A) corrections on the FF distributions from the 
^■^^U (n, /) reaction at i^n = 4 and 5 MeV. The increase 
of D (A) as a function of E^ was under focus and it was 
questioned how u (A) changes. Two assumptions were 
used, the first increasing D (A) for all fragments, the sec- 
ond increasing P (A) for the heavy fragments only. The 
impact found was large on both the yield- and energy dis- 
tributions. The differences are mass-dependent and they 
grew as a function of E^. The post- neutron emission 
mass yields were affected by up to 20-30 % for certain 
masses. Many fission yields in the data libraries were 
measured via the 2E-technique assuming an average in- 
crease of D (A) . If the HE-method proves valid, all these 
yields need to be corrected. Hence, this study emphasizes 
the need of determining 9 (A) as a function of mass and 



V. OUTLOOK 

The important remaining question is how the extra 
neutron emission depends on the fragment mass. The 
Uppsala group and the JRC-IRMM are involved in three 
different projects which may contribute to solving this 



<) Fission yields at IGISOL: Currently, the Uppsala 
group is designing a mono-energetic neutron source 
for the IGISOL- JYFLTRAP facility. The inde- 
pendent fission products will be identified via an 
ion-guide technique and a Penning trap, providing 
a high-precision A/q determination. The aim is 
to study the yield distribution as a function of E^^ 
and to estimate the change of D (A) [T7] . 

<) Coincidence measurements: The JRC-IRMM is 
developing a neutron-detector array to measure 
prompt-fission neutrons in coincidence with FF via 
an ionisation chamber }18 . It is straight forward 
then to measure at different E^ and study how 9 (A) 
changes for the light and heavy fragments. 



The VERDI spectrometer: At the JRC-IRMM, a 
double (v, E) fission-fragment time-of-fiight spec- 
trometer (VERDI) is under development [TH]. By 
measuring both the FF velocities and energies, one 
can estimate D [A) on an event-by-event basis. 
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